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Abstract 

We consider a generalized statistical mechanics model for the cre- 
ation process of cosmic rays which takes into account local tempera- 
ture fluctuations. This model yields Tsallis statistics for the cosmic 
ray spectrum. It predicts an entropic index q given by q = 11/9 at 
largest energies (equivalent to a spectral index of a = 5/2), and an ef- 
fective temperature given by where kTjj 180 MeV is the Hage- 
dorn temperature measured in collider experiments. Our theoretically 
obtained formula is in very good agreement with the experimentally 
measured energy spectrum of primary cosmic rays. 



More general versions of statistical mechanics, as introduced by Tsallis [1 
and further developed by many others [2113], have recently been successfully 
applied to a variety of complex physical systems. The idea is to maximize 
more general entropy measures than the Shannon entropy, which depend on 
a parameter q and which lead to generalized versions of statistical mechanics. 
Ordinary statistical mechanics is contained as a special case for q — 1. In- 
teresting recent physical applications of the so-called nonextensive formalism 
obtained for q ^ 1 include the statistics of fully developed hydrodynamic tur- 
bulence j3], defect turbulence jS], scattering processes in e + e~ annihilation 
OlZj) heavy ion collisions jH], hadron collisions [9,, and models of vacuum 
fluctuations at the Planck scale |lL)j . There is growing evidence that the 
nonextensive formalism, though possessing the mathematical structure of an 
equilibrium formalism, is physically often relevant for nonequilibrium sys- 
tems with a stationary state that possess strong fluctuations of an intensive 
parameter [TT 1 IT2 1 ITS] . 

Very recently Tsallis and Borges [2] analysed the experimentally mea- 
sured energy spectrum of cosmic rays from a nonextensive point of view. 
They showed that a good fit of the measured spectrum can be obtained by 
assuming the interplay of two different nonextensive canonical distributions 
with two different entropic indices q\ and q2 and two different inverse temper- 
atures j3i and j3 2 - Moreover, Kaniadakis jT3] showed that yet another type of 
generalized statistical mechanics (based on entropies that are neither Shan- 
non nor Tsallis) also yields a good fit of the measured spectrum, using three 
fitting parameters. Both these papers represent new interesting ideas, but 
at the same time they use various fitting parameters which are not predicted 
from first principles. The question arises whether one can construct a gener- 
alized statistical mechanics model of cosmic rays that predicts the relevant 
parameters from first principles and that at the same time well reproduces 
the measured cosmic ray energy spectrum using these parameters. 

In this letter we consider such a generalized statistical mechanics model. 
The model is quite generally of relevance for particles that are created by 
scattering processes at very large energies (as cosmic ray particles are). The 
basic idea is that at very high energies the effective thermodynamic inter- 
action volume is small and hence one expects that there are strong local 
temperature fluctuations. A relevant thermodynamic theory of high-energy 
scattering processes is the Hagedorn theory [^Sj, which we will generalize 
here in the sense that we take into account local temperature fluctuations. 
Under reasonable assumptions on the form of the temperature fluctuations 
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this model predicts Tsallis statistics for the cosmic ray spectrum. More 
importantly, it yields a concrete prediction of the entropic index q (or the 
differential spectral index a of the cosmic flux) and of the relevant effective 
temperature that enters into the generalized canonical distributions. Our 
model turns out to yield a very good fit of the experimentally measured cos- 
mic ray spectrum. The advantage of our model is that the most relevant 
parameters are calculated, not fitted. The more subtle features of the cosmic 
ray spectrum (such as the well-known 'knee' and 'ankle') are then finally re- 
lated to the existence of two different populations of cosmic rays, each being 
described by the above nonextensive theory (with the same parameters) but 
having a different flux rate. 

To start with, let us first have a look at the experimental data as collected 
from various experimental groups J7j. These data are displayed in Fig. 1. 
Also shown is a curve that corresponds to a prediction of nonextensive statis- 
tical mechanics. Up to energies of 10 16 eV, the measured flux rate of cosmic 
ray particles with a given energy is well fitted by a generalized canonical 
distribution of the form 

p(E) = C = — . (1) 

(l+/3(g- 

Here E is the energy of the particles 



E = Jc 2 pl + c 2 p 2 y + c 2 p 2 z + m 2 c\ (2) 

(3 = (kT)' 1 is an effective inverse temperature variable, and C is a constant 
representing the total flux rate. For relativistic particles the rest mass m can 
be neglected and one has 

Ep*c\p\. (3) 

The distribution (0) is a g-generalized relativistic Maxwell-Boltzmann distri- 
bution in the formalism of nonextensive statistical mechanics. These kind of 
distributions can be directly obtained by maximizing the Tsallis entropies £Q 

= tttC 1 - £rf) ( 4 ) 

i 

and multiplying with the available phase space volume [SJ E]. The Pi 
are the probabilities of the microstates i. As seen in Fig. 1, the cosmic ray 
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spectrum is very well fitted by the distribution (JTJ) if the entropic index is 
chosen as 

q = 1.215 (5) 
and if the effective temperature parameter is given by 

kf = /T 1 = 107 MeV. (6) 

Our main aim in the following is to derive the form of the distribution (0) from 
physically reasonable assumptions, based on ordinary statistical mechanics 
but leading effectively to a nonextensive description. Moreover, we will derive 
concrete numerical values for the two parameters q and j3, thus proceeding 
from a mere fit towards a theory. 

Cosmic ray particles are created in high energy collision processes as in- 
duced by strong galactic and extragalactic 'accelerators' (see e.g. [THj for a 
review). These 'accelerators' are astrophysical sources and may produce a 
center of mass energy Eqms of the collision process in the GeV and TeV 
region, but much larger energies are possible as well. In general, from a 
statistical mechanics point of view it is well known that the largest aver- 
age temperature that can be reached in high energy collision processes is 
the Hagedorn temperature Th [TO] , experimentally measured to be about 
180 ± 30 MeV [T5J The Hagedorn temperature is much smaller than 
Ecms an d represents a kind of 'boiling temperature' of nuclear matter at the 
confinement phase transition. Even largest Eqms cannot produce a larger 
average temperature than Th due to the fact that the number of possible 
particle states grows exponentially. 

Let us now provide a plausible physical argument how nonextensive be- 
havior comes into play at very large energies. This argument is just based 
on ordinary statistical mechanics but takes into account fluctuations. The 
larger the center of mass energy of the collision process, the smaller the 
volume probed, due to the uncertainty relation. This means, the effective in- 
teraction volume where a thermodynamic description of the collision process 
makes sense will become smaller and smaller with increasing Ecms- How- 
ever, a smaller volume means larger temperature fluctuations. It thus makes 
sense to consider a generalized Hagedorn theory which takes into account 
local temperature fluctuations. 

Assume that locally some value of the fluctuating inverse temperature (3 
is given. We then expect the momentum of a randomly picked particle in 
this region to be distributed according to the relativistic Maxwell-Boltzmann 
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Z{(3) = I E 2 e~ /3E dE = ±. (8) 
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distribution 

vim = -^f^ E - ( 7 ) 

Here p{E\(3) denotes the conditional probability of E given some value of f3. 
We neglect the rest mass m so that E = c\p\. The normalization constant is 
given by 

2 

J 3 

Now let us take into account local temperature fluctuations in the small 
interaction volumes where cosmic ray particles are produced. We have to 
consider some suitable probability density /(/?) of the inverse temperature in 
the various interaction volumes. While the most natural probability distri- 
bution of a continuous random variable taking positive and negative values 
is the Gaussian distribution, the most natural distribution of a continuous 
random variable living on a positive support (such as (3) is the ^-distribution 
of degree n, i.e. the probability density of (3 is given by 

n 

The ^ -distribution is a typical distribution that naturally arises in many cir- 
cumstances, for example if n independent Gaussian random variables Xi, i — 
1, . . . , n with average are squared and added. If we write 



/5:=E X ' ( 10 ) 



i=l 



then p has the probability density function Q. The average of the fluctuating 
j3 is given by 

poo 

(P)=n(X?)= / pf(p)dp = p (11) 



Jo 

and the variance by 

{P 2 )- Pi = 1^1 (12) 

The observed cosmic ray distribution at the earth does not contain any 
information on the local temperature at which the various particles were pro- 
duced. Hence we have to average over all possible fluctuating temperatures, 
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obtaining the measured energy spectrum as the marginal distribution 



p(E) = / p{E\(3)f{(3)d(3. (13) 
Jo 

The integral (fEIJ) with given by Q and p{E\(3) given by (J7J) is easily 
evaluated and one obtains 

p(E) ^ (14) 

(l + (3( q -l)E) — 

with 

q = 1 + (15) 
n + 6 

and 

^ = 4^5- < 16 > 

Eq. (jT3)l is just the fitting function used in Fig. 1. 

The variables Xi describe the independent degrees of freedom contribut- 
ing to the fluctuating temperature. At very large center of mass energies, 
the interaction region is very small, and all relevant degrees of freedom are 
basically represented by the 3 spatial dimensions into which heat can flow. 
We may physically interpret Xf as the heat loss in the spatial i-direction, 
% = x, y, z, during the collision process that generates the cosmic ray particle. 
The more heat is lost, the smaller is the local kT, i.e. the larger is the local 
P given by (fTU|). The 3 spatial degrees of freedom yield n = 3 or, according 

to (USD, 

q = y = 1-222. (17) 

For cosmic rays E C ms is ver Y large, hence we expect a q- value that is close 
to this asymptotic value. The fit in Fig. 1 in fact uses q = 1.215, which 
agrees with the predicted value in eq. (|T7|) to about 3 digits. It also coincides 
well with the fitting value q = 1.225 used by Tsallis et al. Jl] using multi- 
parameter generalizations of nonextensive canonical distributions. 

For smaller center of mass energies, the interaction region will be bigger 
and more effective degrees of freedom within this bigger interaction region 
can contribute to the fluctuating temperature. Hence we expect that for 
finite Ecms n w ih be a bit larger than 3, or q will be slightly smaller than 
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11/9. A good parametrization, in agreement with experimentally measured 
cross sections in e + e~ annihilation, is j7j 

H — e ~ E CMs/E 

^ms) = Q - e . EcMs/Eo (18) 

with E « 45.6 GeV. 

Let us now predict a value for the other 'fitting parameter' entering into 
the nonextensive canonical distribution, the effective temperature kT rs 107 
MeV. The Hagedorn temperature is experimentally measured to be about 
kTu ~ 180 MeV 120] , which is the right order of magnitude. But why is 
the effective temperature fitting the cosmic ray spectrum smaller by a factor 
of about 0.6 as compared to the Hagedorn temperature? 

A possible answer could be related to the way the Hagedorn temperature 
Tjj is experimentally determined. The experimental measurements are eval- 
uated under the assumption that ordinary relativistic Maxwell-Boltzmann 
(MB) statistics applies, which is a good approximation for relatively small 
Ecms- Assuming MB statistics 

p(E) = CE 2 e^ E , (19) 

the temperature (3~ l can be estimated from the maximum E* of the distri- 
bution as 

kT MB = _ E * ( 20 ) 

This relation immediately follows from evaluating p'(E*) = 0. However, 
suppose that in reality at sufficiently large Ecms Tsallis statistics 

p(E) = C (21) 

(l + (q-l)(3E)^ 

applies, due to the fluctuations mentioned before. Then the relation between 
the maximum and the temperature changes. One obtains from evaluating 
p'(E*) = the relation 

p- 1 = kf = ^-y^£* = (3 - 2q)kT MB . (22) 

In collider experiments, the value of the maximum E* is seen to be essentially 
independent of Ecms [Zj • An experimentally determined Hagedorn tempera- 
ture (under the assumption of MB statistics) of 180 ± 30 MeV and a value of 
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q = 1.215 thus yields the effective temperature kf = 0.57-kT MB = (103±17) 
MeV. This just coincides (within the error bars) with the observed effective 
temperature that fits the data in Fig. 1. It also coincides with one of the two 
temperature fitting parameters (96.15 MeV) used by Tsallis et al. in |14j . 
For Ecus oo one obtains from eq. (|T7j) and (J2~2~|) the relation 

kf = -kT 
9 

In this way all relevant parameters determining the shape of the cosmic 
ray energy spectrum are predicted from a simple theoretical consideration 
(except of course for the total flux rate of the cosmic rays, represented by 
the constant C). 

To explain the 'knee' and 'ankle' in the cosmic ray spectrum we adopt 
the standard view [18 . It is reasonable to assume that the knee at E « 10 16 
eV is simply produced by the fact that one has reached the maximum energy 
scale to which typical galactic accelerators can accelerate. This naturally 
implies a rapid fall in the number of observed events with a higher energy, 
i.e. a steeper slope in Fig. 1 between about 10 16 and 10 19 eV. The 'ankle' at 
E 10 19 eV is then produced by the fact that a higher energy population of 
cosmic ray particles takes over from a lower energy population. This higher 
energy population may have a different origin (for example, extragalactic ori- 
gin). The new population has a much smaller flux rate but can reach much 
larger energies. Naturally, the cosmic accelerators underlying the production 
process of this new species of cosmic rays must have a much larger center of 
mass energy Ecms than the ankle energy ~ 10 19 eV, so q should be quite 
precisely given by its asymptotic value 11/9, whereas the effective tempera- 
ture T should be the same as before. The dashed line in Fig. 1 corresponds 
to our formula with q = 11/9, kT = 107 MeV and a flux rate that is smaller 
by a factor 1/50 as compared to the high-flux generation of cosmic rays. It 
is consistent with the data. 

Our nonextensive theory predicts power-law behavior of the measured 
energy spectrum. For large E one has p{E) ~ E~ a where the differential 
spectral index a is given by 

a = J—-2. (24) 
q-1 

q = 1.215 implies a = 2.65 (for moderately large energies), whereas for 
largest energies the asymptotic value q = 11/9 implies a = 5/2. As shown 
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(23) 



in Fig. 2, the largest-energy events are compatible with such an asymptotic 
power law exponent, though of course the error bars of the experimental data 
are huge. 

To conclude, our simple Hagedorn model based on ordinary statistical 
mechanics with local fluctuations of temperature very well reproduces the 
experimentally measured cosmic ray spectrum. By averaging over the tem- 
perature fluctuations one effectively ends up with a nonextensive theory, 
which describes the statistics of cosmic ray particles in a generalized sta- 
tistical mechanics setting. Our model allows for concrete predictions of the 
relevant entropic index q and of the relevant effective temperature parame- 
ter T, and these predicted numerical values are in good agreement with the 
measured data. 
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Fig. 1 Measured energy spectrum of primary cosmic rays (in units of 
m~ 3 s -1 sr^GeV^ 1 ) as listed in ^7j. The solid line is the prediction ((TJ with 
q = 1.215, kT = 107 MeV and C — 5 ■ 10~ 13 in the above units. The dashed 
line is eq. with q = 11/9, kT = 107 MeV and C smaller by a factor 1/50. 
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Fig. 2 Measured cosmic ray energy spectrum E 3 ■ dN/ dE at largest ener- 
gies (data from [THl I2H 122 ) • The straight line is the power law prediction 
with exponent a = 5/2 (corresponding to q — 11/9). 
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